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Abstract 

A simple kinematic model based on superposition p + p collisions, relativistic geometry and 
final-state hadronic rescattering is used to calculate various hadronic observables in ^/sNN = 200 
GeV Au + Au collisions and ^/sNN = 5.5 TeV Pb + Pb collisions. The model calculations are 
compared with experimental results from several ^/sN^ = 200 GeV Au + Au collision studies. If a 
short hadronization time is assumed in the model, it is found that this model describes the trends 
of the observables from these experiments surprisingly well considering the model's simplicity. 
This also gives more credibility to the model predictions presented for ^/sNN = 5.5 TeV Pb + Pb 
collisions. 
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I. INTRODUCTION 



The experiments at the Relativistic Heavy Ion Colhder (RHIC) have produced many 
interesting studies of hadronic observables from relativistic heavy-ion coUisions over the past 
six or so years. The goal has been to use these observables to characterize the conditions 
of the early state of matter in heavy-ion collisions so as to be possible signatures of exotic 
states, such as Quark Matter Q, S S, li I ■ Hadronic observables measured at RHIC can be 
placed into four general categories: spectra, elhptic flow, femtoscopy, and high pt- Examples 
of observables in each category relevant to the present work are the following: "spectra" 



encompasses rapidity, transverse momentum and transverse mass distributions [5|, |6|, 



"elliptic flow" , characterized by the quantity V2, includes V2 vs. r] and V2 vs. px distributions 



0, y, Q: 



11| : "femtoscopy" , also known as Hanbury-Brown-Twiss interferometry 12|. 

3; 
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includes two-pion correlation studies vs. transverse momentum and azimuthal angle 
and "high pt" , which is targeted to be sensitive to jet effects, includes Raa vs. pt and 
dn/dA(f) distributions 0, [3]. Models which describe the early stages of the collision after 
the initial nuclei have passed through each other in terms of partonic degrees of freedom, 
for example as a cascade or in terms of hydrodynamics, have been successful in describing 
the experimental systematics of some of these observables in some kinematic ranges, but 
no single model has thus far succeeded in making an adequate overall description of the 
systematics of all of these observables in a wide kinematical range (for a comprehensive 
overview of comparisons of models with RHIC experiment, see References and 
references therein). 

The goal of the present work is to see how far one can get in describing the experimental 
systematics of all of the observables mentioned above in a wide kinematical range using a 
simple kinematic model with hadronic degrees of freedom. In essence the model is, for each 
heavy-ion collision, a superposition of p + p collisions in the geometry of the colliding nuclei 
with a proper time for hadronization determining the initial space-time position of each 
produced particle, followed by a Monte Carlo hadronic rescattering calculation. The p + p 
collisions are generated by the PYTHIA code [2| at the beam energy of interest. Some of 
the advantages of using this scheme are clear: 

• One has access to all of the particle types available in PYTHIA. 
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• It is conceptually simple: p+p superposition+simple geometry+hadronic rescattering. 

• It has few free parameters - even the hadronization proper time can be set by a 
Tevatron study (see below). 

• Jets are automatically included in the model since PYTHIA has jets and thus hadronic 
observables depending on jets can be studied. 

• The model should be easily scalable via PYTHIA to higher energies such as will 
be found at the Large Hadron Collider (LHC) thus giving the possibility of making 
predictions for heavy-ion collisions at those energies. 

There is no a priori reason why such an approach should be successful, and in fact there 
are reasons to think it should be unsuccessful, the most serious one being that it is hard to 
imagine that hadronic degrees of freedom, rather than partonic degrees of freedom, can be 
valid soon after the nuclei have passed through each other due to the expected high energy 
density. This would require a very short hadronization time in these collisions. On this 
point, it is encouraging that a recent study of pion HBT in Tevatron collisions has shown 
that a similar model for p + p collisions can explain the px and multiplicity dependences for 
the extracted radius parameters if a very short proper time for hadronization of 0.1 fm/c is 
assumed • 

Previous studies using a model similar to this in which hadronic degrees of freedom were 
assumed in the early stage of the heavy-ion collision followed by hadronic rescattering have 
)een shown to give qualitative agreement with experimental results for some observables 
Although similar, there are significant differences between those previous studies 
and the present one: 

• In the old model the initial kinematic state of the hadrons in the collision was pa- 
rameterized as a thermal distribution in px and a Gaussian distribution in rapidity in 
which the temperature of the thermal distribution and width and mid-rapidity den- 
sity of the rapidity distribution (as well as the rapidity densities of different particle 
species) were fixed by comparisons with experiment. In the present model, superposed 
PYTHIA p+p collisions provide all of the information about the initial kinematic state 
of the hadrons, including jets which were not present in the old model. 
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• In the old model the initial hadronic geometry was taken to be similar to a "Bjorken 
tube" , in that there was no initial expansion in the direction transverse to the beam 
direction but initial expansion could occur along the longitudinal direction controlled 
by the hadronization proper time. In the present model, initial expansion is also 
allowed in the transverse direction to satisfy causality in this picture (see Eqs. (1) and 
(2) below). 

• In the previous studies, calculations were done at fixed impact parameter, so com- 
parisons with experimental results which were measured in centrality windows from 
multiplicity cuts were only qualitative. In the present study the model is run in a 
"minimum bias" mode in which a distribution of impact parameters is calculated and 
comparisons with experiments are made using multiplicity cuts to obtain centrality 
windows equivalent to those from experiment, allowing quantitative comparisons with 
experiment. 

• In the previous studies hadronic observables were only calculated for the "soft sector" , 
i.e. pt < 2 GeV/c due to pt distributions for hadrons from the old model becoming 
exponentially larger than experiment for > 2 — 3 GeV/c. It has been found that 
this behavior was due to an error in the inelastic scattering algorithm in the old model 
and has been corrected in the present model. Thus, "hard sector", i.e. > 2 GeV/c, 
studies are now possible. 

Still, the only way to determine whether such a radical and simple picture for heavy-ion 
collisions is valid at all is to compare the results of the model with a range of experimental 
data. Being such a simple model, the main hope will be to give, at best, a qualitative 
description of the trends of the experimental hadronic observables mentioned above. This 
would already be a useful result since it would help establish the hadronization timescale in 
RHIC heavy-ion collisions. 



To this end, model ca^ 



culations will be compared with resn 



ts from the RHIC experiments 



PHOBOS fl, 3, STAR 2, fl ID, Q, 14, 15] and PHENIX [?, Q for Au + Au collisions at 
-y/sjvTv = 200 GeV. The goal will be to make as quantitative comparisons as possible between 
model and experiments. Predictions from the model for LHC-energy Ph + Ph collisions at 
y/SNN = 5.5 TeV will also be given. 
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The paper is organized into the following sections: Section II gives a description of the 
model, Section III presents results of the model for Au + Au collisions at ^/sNN = 200 GeV 
and comparisons with RHIC experiments, Section IV presents predictions from the model 
for ^/sNN = 5.5 TeV Pb + Pb collisions, and Section V gives a summary and conclusions. 



II. DESCRIPTION OF THE MODEL 



The model calculations are carried out in five main steps: A) generate hadrons in p + p 
collisions from PYTHIA, B) superpose p+p collisions in the geometry of the colliding nuclei, 
C) employ a simple space-time geometry picture for the hadronization of the PYTHIA- 
generated hadrons, D) calculate the effects of final-state rescattering among the hadrons, 
and E) calculate the hadronic observables. These steps will now be discussed in more detail. 



A. Generation of the p + p collisions with PYTHIA 

The p + p collisions were modeled with the PYTHIA code lid], version 6.409. The parton 
distribution functions used were the same as used in Ref. [l7j. Events were generated in 
"minimum bias" mode, i.e. setting the low-pr cutoff for parton-parton collisions to zero (or 
in terms of the actual PYTHIA parameter, ckin{3) = 0) and excluding elastic and diffractive 
collisions (PYTHIA parameter msel = 1). Runs were made both with ^/s = 200 GeV and 5.5 
TeV to simulate RHIC and LHC collisions, respectively. Information saved from a PYTHIA 
run for use in the next step of the procedure were the momenta and identities of the "direct" 
(i.e. redundancies removed) hadrons (all charge states) vr, K, p, n, A, p, u, r], r]', 0, and K*. 
These particles were chosen since they are the most common hadrons produced and thus 
should have the greatest effect on the hadronic observables in these calculations. Although 
As were included in the rescattering process (see below) during which they are produced 
abundantly, they were not input directly from PYTHIA since their lifetime is short, i.e. 
^ 1 fm, and thus their initial presence was judged to not effect the results significantly. 
Figure [1] shows an absolute comparison between the invariant cross section distribution for 
charged particles from the PYTHIA p+p run (with resonances decayed) used to generate the 
Au + Au collisions in the present model with a paramerization for 200 GeV p + p collisions 
from PHENIX j^]. As seen, PYTHIA agrees quite well with the PHENIX parameterization 
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FIG. 1: da/dpT from PYTHIA compared with a PHENIX parameterization for 200 GeV p + p 
collisions. 



up to about 6 GeV/c, showing that the p + p collisions input into the model are quite 
reasonable. 



B. Superposition of p + p events to simulate heavy-ion collisions 

An assumption of the model is that an adequate job of describing the heavy-ion collision 
can be obtained by superposing PYTHIA-generated p-\-p coUisions calculated at the beam 
^/s within the collision geometry of the colliding nuclei. Specifically, for a collision of impact 
parameter 6, if f{h) is the fraction of the overlap volume of the participating parts of the 
nuclei such that f{h = 0) = 1 and f{h = 2R) = 0, where R = 1.2A^^^ and A is the mass 
number of the nuclei, then the number of p + p collisions to be superposed will be f{b)A. 
The positions of the superposed p + p pairs are randomly distributed in the overlap volume 
and then projected onto the x — y plane which is transverse to the beam axis defined in the 
z-direction. The coordinates for a particular p + p pair are defined as Xpp, i/pp, and Zpp — 0. 
The positions of the hadrons produced in one of these p+p collisions are defined with respect 
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to the position so obtained of the superposed p + p colhsion (see below). 

For fine tuning of the model so as to get slightly better agreement with the PHOBOS 
rapidity distributions sl , a lower multiplicity cut on the p + p events which are used in 
the heayy-ion calculations was made. For ^s^vat = 200 GeV Au + Au calculations, the cut 
was set to 20, which cut out approximately 26% of p + p events, and for ^/SNN = 5.5 TeV 
Pb + Pb calculations, the cut was set to 38 to cut out a similar fraction of p + p events 
at that energy. It should be emphasized that this was indeed a "fine tuning" cut - if it is 
not made, the overall results presented below from the model do not change appreciably, 
showing the robustness of the model calculations to this cut. 



C. The space-time geometry picture for hadronization 

The space-time geometry picture for hadronization from a superposed p + p collision 
located at {xpp, Upp) consists of the emission of a PYTHIA particle from a thin uniform 
disk of radius 1 fm in the x — y plane followed by its hadronization which occurs in the 
proper time of the particle, r. The space-time coordinates at hadronization in the lab frame 
{xh,yh, Zh,th) for a particle with momentum coordinates {px,Py,Pz), energy E, rest mass 
mo, and transverse disk coordinates {xo,yo), which are chosen randomly on the disk, can 
then be written as 



Xh = Xpp + Xo + T— (1) 

mo 

Vh = Vpp + yo + 'T— (2) 

mo 

Zh = T— (3) 

mo 

E 

th = r— 4 
mo 

Eqs. (1) and (2) show the initial expansion in the transverse direction now present in 
the model. The simplicity of this geometric picture is now clear: it is just an expression 
of causality with the assumption that all particles hadronize with the same proper time, r. 
A similar hadronization picture (with an initial point source) has been applied to — e~ 
collisions |2l|. For all results presented in this work, r will be set to 0.1 fm/c to be consistent 
with the results found in the Tevatron HBT study mentioned earlier 17 1. 



7 



D. Final-state hadronic rescattering 



The hadronic rescattering calculational method used is similar to that employed in previ- 
ous studies Q, Q, except, as mentioned above, the error found in the algorithm to calculate 
inelastic scattering has been corrected. Rescattering is simulated with a semi-classical Monte 
Carlo calculation which assumes strong binary collisions between hadrons. Relativistic kine- 
matics is used throughout. The hadrons considered in the calculation are the most common 
ones: pions, kaons, nucleons and lambdas (vr, K, N, and A), and the p, u;, r], 77', 0, A, and K* 
resonances. For simplicity, the calculation is isospin averaged (e.g. no distinction is made 
among a vr^, vr^, and vr^). 

The rescattering calculation finishes with the freeze out and decay of all particles. Starting 
from the initial stage (t = fm/c), the positions of all particles in each event are allowed to 
evolve in time in small time steps (At = 0.5 fm/c) according to their initial momenta. At 
each time step each particle is checked to see a) if it has hadronized {t > t^, where th is given 
in Eq. (4)), b) if it decays, and c) if it is sufficiently close to another particle to scatter with 
it. Isospin-averaged s-wave and p-wave cross sections for meson scattering are obtained from 



Prakash et al. 



22I ] and other cross sections are estimated from fits to hadron scattering data 



in the Review of Particle Physics 23|] . Both elastic and inelastic collisions are included. The 
calculation is carried out to 200 fm/c for a/J/vtv = 200 GeV Au + Au collisions and to 400 
fm/c for y/SNN = 5.5 TeV Pb + Pb collisions which allows enough time for the rescattering 
to finish (as a test, calculations were also carried out for longer times with no changes in the 
results). Note that when this cutoff time is reached, all un-decayed resonances are allowed to 
decay with their natural lifetimes and their projected decay positions and times are recorded. 

Figure [2] shows the time evolution plotted up to 50 fm/c of the particle density calculated 
at mid-rapidity, i.e. in the rapidity range —l<y<l, and the number of rescatterings 
per time step from the model for minimum bias (see below) ^/snn = 200 GeV Au + Au 
collisions. The solid lines show the average values of the quantities whereas the dotted lines 
show the average +a to give an idea of the width of the distribution. The density is seen to 
start out high at 6-10 fm~'^ and then to fall off rapidly with time such that by 4 fm/c the 
density is at or below 1 fm^^. The rescattering rate starts small in the first time bin due 
to time dilation of the hadronization time and the requirement imposed in the model that 
particles must hadronize before they can scatter. By the second time step the rescattering 
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FIG. 2: Time evolution up to 50 fm/c of the particle density calculated at mid-rapidity (—1 < y < 
1) and the number of rescatterings per time step from the model for minimum bias y/s^N = 200 
GeV Au + Au collisions. The solid lines show the average values of the quantities whereas the 
dotted lines show the average +a. 



rate increases quickly and then falls off rapidly with time as does the density. The time 
evolution of the density will be discussed more later. 

The rescattering calculation is described in more detail elsewhere 
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19|. The validity 



of the numerical methods used in the rescattering code have recently been studied using the 



subdivision method, the results of which have verified that the methods used are valid 



20|. 



E. Calculation of the hadronic observables 



Model runs are made to be "minimum bias" by having the impact parameters of collisions 
follow the distribution da/db oc b, where < b < 2R. Observables are then calculated from 
the model in the appropriate centrality bin by making multiplicity cuts as done in the 
experiments, as well as kinematic cuts on rapidity and pt- For the present study, a single 
87K event minimum bias run was made from the model for y/SNN = 200 GeV Au + Au 
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collisions which was then used to calculate all of the hadronic observables from the model 
which are shown in this work for that system. In this way, a consistent picture of the 
agreement between the model and experiments can emerge since it is virtually impossible to 
optimize the model to agree with experiment for a particular observable without spoiling the 
agreement for others which are calculated from the same run. In the case of the y^J/viv = 5.5 
TeV Pb + Pb predictions, a single minimum bias run with 800 events was made. 



III. RESULTS FROM THE MODEL FOR Au + Au COLLISIONS AT ./sWn = 200 
GEV AND COMPARISONS WITH RHIC EXPERIMENTS 

Various hadronic observables have been calculated from the 87K minimum bias run from 
the model mentioned above and are now compared with measurements from RHIC experi- 
ments. The observables and experiments with which to be compared are: 

• Spectra - dn/d?] (PHOBOS), dn/dpr (PHENIX), dn/dniT (STAR) 

• Elliptic flow - V2 vs. Tj (PHOBOS), V2 vs. pr charged particles (STAR), V2 vs. pr 
identified particles and Vijuq vs. px/nq (PHENIX) 

• HBT - TTTT VS. and vrvr vs. kx (STAR) 

• High pt - Raa vs. Pt (PHENIX), dn/dA(f) vs. (STAR) 

In the spirit of making as quantitative comparisons as possible between model and experi- 
ments, unless explicitly specified otherwise, absolute normalizations are used for the model 
observables in the plots shown. 



A. Spectra 

Figures [31 12 and O show model comparisons with PHOBOS , PHENIX 0, and STAR 
0] for 7], Pt, and tjit distributions, respectively. In Figures [3] and H] the centrality dependence 
of charged hadrons is also shown, and in Figure [5] the particle species dependence is shown. 
In Figure [3], the centrality dependence of the rapidity distribution is followed fairly well by 
the model, although it is seen that the model shapes are slightly broader than experiment 
and the overall agreement near mid-rapidity is at the 10-15% level. Since the model is 
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"isospin averaged", the model distributions are multiplied by 2/3 to approximate all charged 
particles. In Figure HI the experimental centrality dependence is once again seen to be 
described reasonably well by the model, especially the absolute scale at low px- At higher 
Pt, particularly in the minimum bias case, the model under predicts the absolute scale by 
as much as a factor of four at some points over the pt range, yet still follows the trend 
of the data up to the highest pt shown. To approximate {h^ + h~)/2 for the model, the 
model distributions are multiplied by 1/3. In Figure[5]the absolute scale of the experimental 
particle species dependence is reproduced reasonably well by the model, as well as the "radial 
flow effect" of the slopes decreasing with increasing particle mass, but the absolute slopes 
from the model at low rriT — ttlq are seen to be uniformly somewhat too large (as is also 
seen in Figure H] at low px)- The model distributions are multiplied by 1/3 to approximate 
positive charges. 

Although there are clear differences seen in details between model and experiments as 
described above, it is still remarkable that this simple model does as well as it does in repro- 
ducing the overall absolute scales and dependences of these "bread and butter" experimental 
distributions. It is judged that the description by the model of these "basic observables" 
is adequate enough to cautiously proceed with using the model to calculate the "derived 
observables" such as elliptic flow, HBT, etc.. 

B. Elliptic flow 

The elliptic flow variable, V2, is defined as 

V2 =< cos(20) > (5) 
(J) = arctanf— ) 

Px 

where "<>" implies a sum over particles in an event and a sum over events and where 
Px and Py are the x and y components of the particle momentum, and x is in the impact 
parameter direction, i.e. reaction plane direction, and y is in the direction perpendicular to 
the reaction plane. The V2 variable is calculated from the model using Eq. ([5]) and taking 
the reaction plane to be the model x — z plane. 

Figures IMTT] show comparisons between the model and experiments for elliptic flow. 
Figure El compares the model to PHOBOS for charged particles for V2 vs. rj in a. centrality 
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FIG. 3: Rapidity distributions for Model compared with PHOBOS for several centralilties. 



window of 25 — 50% [8]. The model is seen to agree with the measurements within error bars 
for the entire range in rj, i.e. —6 < rj < 6, although it looks systematically about 10% lower 
than experiment around mid-rapidity. Note that in the model, V2 is completely determined 
by rescattering such that if the rescattering is turned off, V2 = in all kinematic regions. 



Figure [7| shows a comparison of the mode 
particles in several centrality bins from STAR 



to V2 vs. pt for pt < 2 GeV/c for charged 
9I. The model is seen to do a reasonable job in 



representing the different centralities, although it is systematically higher than measurement 
by about 0.01 — 0.02 over the entire pt range for the — 5% centrality bin. Figure [H] compares 
the model with STAR for V2 vs. pt for charged particles in a centrality bin 10 — 40% in 
a wide pt range, i.e. px < Q GeV/c [10!]. What is remarkable about this comparison is 
that the model describes the pt behavior of the experiment in which V2 increases for p^ < 2 
GeV/c, flattens out, and then starts decreasing forp^ > 3 GeV/c. Once again, this behavior 
is completely rescattering-driven in the model. 

In Figures [9] and [10] the model is compared with identified-particle V2 vs. pt plots for 
caons, and pions and protons, respectively, and for minimum bias centrality from PHENIX 



ll| . For both figures the low px behavior, i.e. px < 2 GeV/c is described well and quan- 
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FIG. 4: pt distribution for charged hadrons for Model compared with PHENIX for several cen- 
tralities. The meaning of "mid-r/" in this case is that the PHENIX spectra are measured in the 
range —0.18 < rj < 0.18, whereas the model spectra are calculated in the range —1 < r] < 1 for 
better statistics. 



titatively by the model, whereas the high pt behavior, i.e. pt > 2 GeV/c is only described 
qualitatively. In Figure [9] the high px behavior of the kaons flattens out for both experiment 
and model, but for the model it flattens out to a slightly lower value. In Figure [10] the high 
Pt behavior for the pions is to flatten out and start decreasing as the model also does but to 
a lower value whereas the proton V2 continues increasing, as it does for the model nucleons, 
but the model does not increase as fast (and the model 3.5 GeV/c point decreases). Figure 
[TT] shows the plots in Figures [9] and [10] replotted in terms of the number of valence quarks 
in the identified particle, Ug, as V2/nq vs. pr/nq. The point of doing this is to show that the 
different identified particles follow a universal curve when plotted on the same graph this 
way. Not surprisingly in the context of the discussion above, the model is seen to follow the 
experimental scaling quantitatively for pr/nq < 1 GeV/c and qualitatively at a lower value 
for Pr/nq > 1 GeV/c. 
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FIG. 5: niT distributions for identified particles for Model compared with STAR. 



C. Two-pion femtoscopy (Hanbury-Brown-Twiss interferometry) 

For the HBT calculations from the model, the three-dimesional two-pion correlation func- 
tion is formed and a Gaussian function in momentum difference variables is fitted to it to 
extract the pion source parameters. Boson statistics are introduced after the rescattering has 
finished (i.e. when all particles have "frozen out") using the standard method of pair-wise 
symmetrization of bosons in a plane- wave approximation [24] . The three-dimensional corre- 
lation function, C{Qside, Qout, Qiong), is then calculated in terms of the momentum-difference 
variables Qside, which points in the direction of the sum of the two pion momenta in the 
transverse plane, Qout, which points perpendicular to Qside in the transverse plane and the 
longitudinal variable along the beam direction Qiong- 

The final step in the calculation is extracting fit parameters by fitting a Gaussian param- 
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eterization to the model-generated two-pion correlation function given by, 

C (^Q side, Qout, Qiong) 

G[l + A exp(— Qsjcie-^side ~ QoutR-out ~ QlonqRlonq ~ Q outQ sideRoutside)] (6) 
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FIG. 6: V2 vs. T] for Model compared with PHOBOS 25-50% centrality. The inner error bars on 
the PHOBOS data points show the statistical error and the outer error bars show the statisti- 
cal+systematic error. 



where the i?-parameters, called the radius parameters, are associated with each momentum- 
difference variable direction, G is a normalization constant, and A is the usual empirical 
parameter added to help in the fitting of Eq. to the actual correlation function (A = 1 
in the ideal case). The fit is carried out in the conventional LCMS frame (longitudinally 
comoving system) in which the longitudinal boson pair momentum vanishes 25| . Figure [T2] 
shows a sample projected correlation function from the model with projected fit to Eq. (I6l). 

Figures [13] and [TH show comparisons between the model and STAR for radius parameters 
extracted as a function of azimuthal angle, 0, for two centrality cuts, — 5% and 40 — 80%, 
respectively |lii]. In FigurefTSlthe model is seen to describe the more or less "flat" dependence 
on for Rout, Rside, and Riong seen in the experiment for these central collisions, although 
the model is seen to under predict the magnitude of Rside by about 20%. The model is 
seen to follow the trend of the experiment for Rlutside within the large statistical error bars 
shown. For the less central collisions shown in Figure [TU the model describes the oscillatory 
behavior now seen in Rout and Rside as well as R^^tside the continued fiat dependence 
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FIG. 7: V2 vs. pt for Model compared with STAR for charged particles and several centralities. 
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FIG. 9: V2 vs. pt for Model vs. PHENIX for kaons. The meaning of "mid r/" in this case is that 
the PHENIX spectra are measured in the range —0.35 < rj < 0.35, whereas the model spectra are 
calculated in the range — 1.3 < r/ < 1.3 for better statistics. 



seen in Riong, although once again under predicting the magnitude of Rside by about 30%. 
The A-parameters extracted in the fits from the model were constant in (p with the values 
0.61 and 0.54 for the — 5% and 40 — 80% centrality bins, respectively. 

Figure [T5l compares the model with STAR for the fc-r dependence of the radius parameters 
in a centrality bin of — 5% [l4|. For these fits, the parameter Routside in Eq. ([6]) is set 
to zero. Several methods used by STAR in extracting their fit parameters were used and 
as seen they give approximately the same basic results jl4]. The model is seen to follow 
the general trend of the experiment for decreasing radius parameters for increasing kr- The 
model slightly over predicts this effect for Riong and, as was the case in the azimuthal HBT 
results shown earlier, consistently under predicts the magnitude of Rside by about 24%. As 
seen in the experiment, the model A-parameter is less than 1, being constant in kx at about 
0.61, but larger than the experimental values on average by about 30%. Note that the main 
source of A < 1 in the model is from the presence of long-lived resonances such as r] and 77'. 
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FIG. 10: V2 vs. pt for Model vs. PHENIX for pions and protons. The meaning of "mid 77" is the 
same as for Figure [H 



D. High Pt 

Studying the high px behavior of the observables Raa and dn/dAcj) is thought to be a way 
of more directly studying QCD processes, such as jets, in heavy-ion collisions. As mentioned 
earlier, since the present model is based on using PYTHIA which uses QCD processes in 
calculating p+p collisions, the model should contain these effects and thus should be suitable 
for comparing with experiments which measure these observables. 

The Raa is defined as 

= ^A^^MT^^^^^"^"^^ 
where the numerator is the usual Au + Au px distribution as shown in Figure H] and the 
denominator has the p + p pr distribution normalized to cross section as in Figure [1] and 
multiplied by the quantity TauAu-, the Glauber nuclear overlap function which is different 
for each centrality cut. 

Figure [16] compares the model to PHENIX for Raa vs. pt for three centrality windows 
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FIG. 11: V2/nq vs. px/riq for Model vs. PHENIX for pions, kaons, and protons. 



0]. For both the model and PHENIX the plots in Figures [T] and H] were used in Eq. ([7]) 
to calculate these Raa plots along with the Taa values shown in Table I. of Ref. [3]. The 
error bars shown for the PHENIX plots are a sum of both statistical error and the overall 
scale uncertainty, and they mostly reflect the scale uncertainty. As seen, the model describes 
three main qualitative features of the experiment: 1) for large pt the Raa decreases with 
increasing pt-, and as the centrality window goes from minimum bias (0 — 92%) to peripheral 
(80 — 92%) 2) the scale of Raa increases and 3) the dependence of Raa on Vt tends to flatten 
out. It is also seen that, even with the uncertainty in the PHENIX overall normalization, 
the model scale tends to be lower than experiment, and at low p^ the peaks in the plots 
for the model occurs at about 1.3 GeV/c whereas the peaks occur at about 2.3 GeV/c for 
experiment. These differences of the model with experiment reflect the differences already 
seen in the pt distributions in Figure HI but the qualitative similarities as discussed above 
are clearly present out to the highest px shown. 

Figure [T7l shows dn/d/S.(j) vs. plots from the model and a comparison of one of them 
with STAR charged particles [isl. The model plots, which include all hadrons, are made 
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FIG. 12: Sample two-pion correlation function with Gaussian fit projected onto the Qout^ Qside, and 
Qiong axes from the Model. The collision centrality is — 5% with cuts on the pions —0.5 < y < 0.5, 
0.15 <pt < 0.8 GeV/c, and 0.15 < kr < 0.25 GeV/c. 

using the same cuts on rapidity and pt as used by STAR, namely for individual particles 
I |< 0.7 and pt > 2 GeV/c, and for particle pairs, one of which is a "trigger particle", 
from which A(j) is formed, | Ai] \< 1.4 and p^"^ > 4 GeV/c. The lines are fits to the model 
points to guide the eye and the model dn/dAcf) normalizations are in counts per bin. Figure 
[T7k shows a plot from PYTHIA for 200 GeV p+p as a reference. The forward and backward 
peaks from di-jet production are clearly seen at A0 = and ivr, respectively. Figure [TTb 
compares the model to STAR minimum bias ^/snn = 200 GeV Au + Au collisions. The 
scale for the STAR plot is shown on the right-hand axis of the figure. As seen the model 
describes the shape of the experiment well, both the width of the forward peak and the 
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FIG. 13: Azimuthal two-pion HBT parameters vs. from Model vs. STAR for centrality 0-5%. 
Pions are accepted in the cut ranges — 1 < y < 1, 0.1 < pr < 0.6 GeV/c, and 0.15 < kx < 0.6 
GeV/c. 



relative height of the backward peak with respect to the height of the forward peak. Figure 
[T7b shows a model plot for a medium centrality window at 10 — 30%. Features similar to the 
minimum bias plot are seen with a forward and backward peak of similar relative heights, 
although the width of the forward peak is somewhat larger than in the minimum bias case. 
A more central case from the model is shown in Figure [TTH . where the centrality window is 
— 10%. Although this plot is pushing the edge of the statistics possible from the 87K event 
model run used in this study for -.Jsnn = 200 GeV Au + Au collisions, it appears to have 
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FIG. 14: Azimuthal two-pion HBT parameters vs. (j) from Model vs. STAR for centrality 40-80%. 
Pions are accepted in the cut ranges — 1 < y < 1, 0.1 < px < 0.6 GeV/c, and 0.15 < kx < 0.6 
GeV/c. 



a qualitatively different shape compared with the other plots shown in this figure. Namely, 
besides the presence of the forward peak, the plot looks more or less ffat for values of A(j) 
out to ±7r, i.e. the backward peak appears suppressed. This is the same general behavior 



seen in STAR in the same centrality window, i.e. Figure Ic of Ref. 
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FIG. 15: Two-pion HBT parameters vs. kx from Model vs. STAR for centrality 0-5%. Pions are 
accepted in the cut ranges —0.5 < y < 0.5 and 0.15 < pr < 0.8 GeV/c. The meaning of the three 
methods used by STAR to extract parameters may be found in Ref. \l£ 



IV. PREDICTIONS FROM THE MODEL FOR Pb + Pb COLLISIONS AT ^/i]^ = 
5.5 TEV 

Predictions from the model for LHC-energy collisions for y/s^N = 5.5 TeV Pb + Pb 
collisions are presented below. These will be sample predictions for px and rj distributions, 
elliptic flow and HBT to give a flavor of the differences predicted by the model for LHC 
Pb + Pb collisions compared with RHIC Au + Au collisions. In this spirit, RHIC-energy 
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FIG. 16: Raa for Model compared with PHENIX for several centralities. 



Au + Au collisions with the same kinematic conditions as the LHC-energy Pb + Pb collisions 
will be shown to indicate the trends of the predictions. With the changes described above 
in Section II, the model will be used in the same way to make the y/sm^ = 5.5 TeV Pb + Pb 
predictions as it was used for the ^/snn = 200 GeV Au + Au calculations, including the use 
of the short proper hadronization time of r = 0.1 fm/c. 

Figures [18] - [21] show the model predictions for ^/sNlv = 5.5 TeV Pb + Pb collisions 
compared with ^/snn = 200 GeV Au + Au collisions, also from the model. The most 
noticeable features of these predictions are summarized below: 

• dn/ dr] near mid-rapidity for charged particles is predicted to be about 1400 for a — 5% 
centrality window at the LHC. - This is seen in Figure [TH which predicts that the 
mid-rapidity charged particle density for LHC Pb + Pb will be about a factor of 2.5 
greater than for RHIC Au + Au. From Figure [3] it was seen that the model density 
was about 10% lower than experiment. Thus even boosting the LHC prediction up 
by 10% puts its value at the lower end of the range of predictions which have been 
recently made of 1500 — 4000 in central collisions using various extrapolations of RHIC 



experimental rapidity densities |26|. 



The charged particle pt distribution for pt > 5 GeV/c is predicted to be about two 
orders of magnitude larger at the LHC compared with RHIC. - This is seen in Figure 
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FIG. 17: dn/dA(p vs. Ac/) plots using STAR cuts on px- a) 200 GeV p+p collisions from PYTHIA. 
b) \/snn = 200 GeV Au + An with minimum bias centrality for Model and STAR, c) ^smn = 200 
GeV Au + Au for Model with 10 — 30% centrality. d) Same as c) with — 10% centrality. Fits to 
the points are also shown to guide the eye. 



[m and is an expected consequence of the higher ^/snn in the LHC collisions that the 
Pt distributions at high px should be greatly enhanced. 

• Elliptic flow in minimum bias centrality collisions is predicted to be slightly smaller 
at the LHC compared with RHIC. - As seen in Figure [20] the plot of V2 vs. pt for all 
hadrons for LHC Ph + Ph looks similar to that for RHIC Au + Au collisions, but overall 
the LHC plot gives a slightly smaller V2 for the entire range in px- This is somewhat 
unexpected since in this model V2 is produced exclusively by the rescattering process, 
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FIG. 18: Rapidity distributions from Model comparing Pb + Pb collisions at ^/sNN = 5.5 TeV 
(LHC) with Au + Au collisions at ^/sNN = 200 GeV (RHIC) for charged particles and — 5% 
centrality. 

and since the rapidity density seen in Figure [18] is higher at LHC, it might be expected 
that the more rescattering would resuh in a higher V2 than for RHIC. Rather than 
this, it seems as though the initial "almond-shaped" geometry present in the 6 > 
collisions starts to become "washed out" from the enhanced number of rescatterings. 
This will be discussed more below. 

• TTTT HBT radius parameters in — 5% centrality collisions are predicted to be 20 — 
30% larger at the LHC compared with RHIC. - In Figure [2T] the behavior of the kT 
dependence of the HBT parameters for LHC Pb + Pb is similar to what is seen for 
RHIC Au + Au, but the overall scale for the radius parameters is slightly larger. 
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FIG. 19: pt distributions from Model comparing LHC Pb + Pb and RHIC Au + Au collisions. 



V. DISCUSSION 

Although some discussion of the individual results has been given above, it is useful to 
reflect on the sources of the effects in the present model which produce these results and 
to describe what interconnections may exist among them. Whereas PYTHIA provides the 
baseline p + p kinematics, the hadronic rescattering is responsible for all collective effects in 
the model beyond p + p. As seen in Figure [2], the rescattering rate is large in the early stage 
of the collision with some rescattering persisting to times of 50 fm/c and beyond. The radial 
flow effects seen in the identified particle distributions of Figure [5] and kx dependence 
of the HBT parameters in Figure [15] as well as the elliptic flow effects seen in the V2 plots of 



Figures [& 



2 are mostly established in the early stage, i.e. by t ~ 10 fm/c (see Figure 18 of 



Reference [18|] for a calculation of y/s^N = 200 GeV Au+Au collisions using a model similar 
to the present one). Comparison of the model with the experimental azimuthal HBT results 
shown in Figures [13] and [H] provides a double test of the flow generated by the rescattering 
since both radial flow effects, i.e. the sizes of the HBT parameters, and elliptic flow effects, 
i.e. their dependence, are present. While the present model gives a reasonably good 
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FIG. 20: V2 vs. pT from Model comparing LHC Pb + Pb with RHIC Au + Au collisions; minimum 
bias centrality, all hadrons, and —1 < r] < 1. 



description of the experimental elliptic flow results using rescattering effects alone, other 
cascade studies have found that either it is not possible to generate large enough elliptic 
flow 271, l28|] or it is necessary to use extreme elastic parton cross sections [29] to describe 
experiments. Although there are many differences in detail between the present study and 
those studies, the basic feature allowing the present model to generate enough elliptic flow 
is that hadronic rescattering is allowed to take place from the earliest times and thus during 
the highest densities, as seen in Figure O A concern regarding carrying out transport 
calculations at high densities is that non-physical superluminal artifacts can be introduced 
which can effect the results [29]. As mentioned earlier, a study of this effect has been carried 
out for a model very similar to the present one and it was found that the observables from 
the model studied, i.e. spectra, elliptic flow, and HBT, were not signiflcantly affected |20 |. 
Thus it is considered unlikely that these artifacts play a signiflcant role in the present study. 

This last statement should also be true for the ^/sNN = 5.5 Tev Pb + Pb predictions 
shown in Figures fT8ll2Tl Although the rapidity density at mid-rapidity for LHC Pb + Pb is 
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FIG. 21: Two-pion HBT vs. kr from Model comparing LHC Pb+Pb with RHIC Au+Au collisions; 
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seen to be more than twice as large as for RHIC Au + An in Figure [181 the mid-rapidity 
particle density for LHC Ph + Ph for early times is found to be similar to that for RHIC 
Au + Au seen in Figure |2l This is due to the larger hadronization volume and time resulting 
from the higher average particle momenta generated in the LHC-energy collisions, as seen in 
Figure fT9l and calculated from Eqs. (l)-(4). This similarity between the particle densities 
would explain why the flow effects seen in Figures [20] and [5T] are similar, i.e the similarity 
between the LHC and RHIC elliptic flow and the dependence of the HBT radius parameters 
on (although the overall size of the LHC radius parameters is greater). 

The features of the high px observables Raa and dn/d^cf) calculated from the model 
and shown in Figures [16] and [TT] are also driven by the underlying rescattering. For these 
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observables it is the energy loss of the p + p produced high pt particles rescattering with 
the rest of the particles that creates the effects. The effect of the rescattering can be seen 
for Raa by considering that if the rescattering were turned off in the model the numerator 
in Eq. [3, which would be simply a superposition of p + p events, would be exactly the same 
distribution as the denominator except for a scale factor, resulting in a flat dependence on 
Pt for all centralities. The same would be the case for dn/dAip in Figure [T71 if rescattering 
were turned off: Figure [TTB . would look like Figure [T7k . 

VI. SUMMARY AND CONCLUSIONS 

A simple kinematic model based on superposition oip + p collisions, relativistic geometry 
and final-state hadronic rescattering has been used to calculate various hadronic observables 
in y/SNN = 200 GeV Au + Au collisions and ^Jsnn = 5.5 TeV Pb + Pb collisions. The 
model calculations were compared with experimental results from several ^^snn = 200 GeV 
Au + Au collision studies from RHIC. With the short hadronization time assumed in the 
model of r = 0.1 fm/c, it is found that this model describes the trends of the observables 
from these experiments surprisingly well considering the model's simplicity. This also gives 
more credibility to the model predictions presented for LHC-energy ^/snn = 5.5 TeV Pb+Pb 
collisions. 

As shown above, the main strength of the present model is not that it gives precise 
agreement with experiment for individual observables in particular kinematic regions, but 
in its ability to give an overall qualitative description of a range of observables in a wide 
kinematic region, i.e. to summarize the gross features seen in experiments for ^/snn = 200 
GeV Au + Au collisions. Another strength is its simplicity. Besides the kinematics generated 
in the superposed p + p collisions by PYTHIA, the only other "active ingredient" in the 
model driving the kinematics underlying the hadronic observables shown is the final-state 
hadronic rescattering. As discussed above, if the hadronic rescattering were turned off in 
the model, all elliptic flow would disappear, the HBT radius parameters would lose all and 

dependence and would be significantly smaller, the Raa vs. pt plots would be fiat, and 
all of the dn/dA(f) plots would look like p + p (i.e. Figure [T7k). Making r large is another 
way to effectively turn off all of these effects since the rescattering is very sensitive to this 
variable since it controls the initial particle density, e.g. for r = 1 fm/c all of these effects 
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would already be greatly reduced [l7j. The price to be paid for this simplicity is to assume 
that either hadrons or "hadron-like" objects can exist in the earliest stage of the heavy- 
ion collision just after the two nuclei pass through each other, i.e. that the hadronization 
time in the frame of the particle is short and insensitive to the environment in which it 
finds itself. Clearly this simple picture is an oversimplification as is seen in the model's 
shortcomings from the comparisons given above with RHIC experiments. However, some 
of these shortcomings could possibly be improved by "perturbatively" adding a few extra 
pieces of physics to this model, but thereby making it less simple. 
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